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ABSTRACT 



Aims. We propose a new, more realistic description of the perturbed gravitational potential of spiral galaxies, with spiral arms having 
Gaussian-shaped groove profiles. The aim is to reach a self-consistent description of the spiral structure, that is, one in which an initial 
potential perturbation generates, by means of the stellar orbits, spiral arms with a profile similar to that of the imposed perturbation. 
Self-consistency is a condition for having long-lived structures. 

Methods. Using the new perturbed potential we investigate the stable stellar orbits in galactic disks for galaxies with no bar or with 
only a weak bar. The model is applied to our Galaxy by making use of the axisymmetric component of the potential computed from 
the Galactic rotation curve, in addition to other input parameters similar to those of our Galaxy. The influence of the bulge mass on 
the stellar orbits in the inner regions of a disk is also investigated. 

Results. The new description offers the advantage of easy control of the parameters of the Gaussian profile of its potential. We 
compute the density contrast between arm and inter-arm regions. We find a range of values for the perturbation amplitude from 400 to 
800 km 2 s~ 2 kpc~', which implies an approximate maximum ratio of the tangential force to the axisymmetric force between 3% and 
6%. Good self-consistency of arm shapes is obtained between the Inner Lindblad resonance (ILR) and the 4:1 resonance. Near the 
4: 1 resonance the response density starts to deviate from the imposed logarithmic spiral form. This creates bifurcations that appear 
as short arms. Therefore the deviation from a perfect logarithmic spiral in galaxies can be understood as a natural effect of the 4:1 
resonance. Beyond the 4: 1 resonance we find closed orbits that have similarities with the arms observed in our Galaxy. In regions near 
the center, elongated stellar orbits appear naturally, in the presence of a massive bulge, without imposing any bar-shaped potential, 
but only extending the spiral perturbation a little inward of the ILR. This suggests that a bar is formed with a half-size ~ 3 kpc by a 
mechanism similar to that of the spiral arms. 

Conclusions. The potential energy perturbation that we adopted represents an important step in the direction of self-consistency, 
compared to previous sine function descriptions of the potential. In addition, our model produces a realistic description of the spiral 
structure, which is able to explain several details that were not yet understood. 

Key words. Galaxies: spiral - Galaxies: structure - Galaxy: kinematics and dynamics - Galaxy: disk - galaxies: bulges 



1. Introduction 

D espite its historica l merit, the spiral structure theory, proposed 
bv lLin&ShvJ(ll964h . has well-known limitations. For instance, 
Lin & Shu add gas and star into a single component, when their 
dynamics differ in response to the same potential; stars, for ex- 
ample, are not affected by pressure gradients. Additionally, real 
arms are not tightly wound, as assumed in that theory, and the 
physical properties of the disk (densities, velocity perturbations) 
do not have a sine-shaped dependence in azimuth , as arg ued in 
the present paper. Around the same time, iToomrel (ll964l) . in an 
opposing point of view, argued that since the largest fraction of 
matter in the disk is in the form of stars, one can, in a first ap- 
proximation, neglect the gas and work with stellar dynamics. He 
devoted much of his attention to understanding why an excess 
of stars in a given region does not grow without limit, due to 
the excess of gravitational attraction. Regardless, his model for 
the arms focused more on stellar dynamics, considering them as 
regions with stellar excess of stars, that rotate with velocity that 
does not diverge from the velocity of the disk matter. 

iKalnais! (1 1972b proposed an innovative way of understanding 
the spiral structure, which was focused solely on stellar dynam- 
ics. He showed that by introducing some degree of organization 



to the closed stellar orbits, we can generate regions where these 
orbits are crowded. Such an organization is achieved by rotating 
orbits of increasing size by a given angle, relative to previous 
ones. These regions of high stellar densities look like logarithmic 
spirals and rotate with the required pattern speed to transform 
the usually open orbits into closed ones. Furthermore, such an 
organ ization of stellar orbits can be achieved via galactic colli- 
sions (lGerber&Larnbll994 . The great advantage of the Kalnajs 
model in comparison to Toomre's is that the former's spiral 
structure is not disrupted by differential rotation; in other words, 
it can be long lived. It is important to note that the search for sta- 
ble or quasi-stationary solutions has traditionally been the aim 
of dynamical models of disks. This is still a legitimate concern, 
given that the near-infrared morphology of several grand design 
galaxies reveals the existence of an u nderlying two-arm struc- 
ture constituted mainly of old stars (|Gr osb0l & Patsis 1998; 
i Block & Wainscoatl 1 199 it iRix & Rieke 1 119931: iRix &Zar itskv 



1 994: IBlock et all 19941) . Also. lLepine et alJ(l201 lbl) argued that, 
based on the metallicity gradient, the present spiral structure of 
Galaxy has an age of at least several billion years. 

The Kalnajs model is the most promising one, as it gives 
a correct understanding of the nature of the arms as regions 
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where stellar orbits are crowded. However, Kalnajs' model is 
only a first approximation, since it makes use of stellar orbits 
associated with an axisymmetric potential. In order to reach 
self-consistency, one must introduce the effect of the perturba- 
tions of the potential on the stellar orbits. The perturbations 
produce changes in the shapes of the orbits, which in turn 
change the regions where the orbits are crowded and gener- 
ate a new perturbation. If at the end of this cycle we obtain 
approximately the same density distribution as at the end of 
the starting one, we say that the model is self-consistent, and 
that the structure will be long lived. Self-consistent models, or 
at least approximately self-co nsistent models in terms of arm 
shapes, were constructed bv IContopoulos & Grosb0l] dl986), 
IContopoulos & Grosbcii (119881) andlPatsis et all d 199 lb for large 

cm i-'i 1 ii'i I 'iv i ai' onn hi; A m orol xt T om r\c± i 1 QQh\ morooftor A xt T \ 



spiral galaxies, and by Amaral & Lepine ( 1996) (hereafter A&L) 
and iPichardo et al.l (120031) for the special case of our Galaxy. In 
the works of Contopoulos and collaborators and of A&L, the 
classical expression for the perturbation Oi was used in a frame 
of reference rotating with constant velocity Cl p : 



<t>i = £oRe~ s ' cos m 



InR 
tan i 



imp 



(1) 



with some diversity in the amplitude £o, the behavior of the am- 
plitude with radius (here represented by R e~ e ' R ), the number of 
arms (m), and the pitch angle i. A property of this expression 
is that if we fix the radius R, the variation of the potential as a 
fu nction of the azimuth al angle is a sine function. In the work 
of Pichaix fo et al.l (120031) . the perturbing potential is represented 
as a superposition of a large number of oblate spheroids, with 
their centers distributed along the locii of the spiral arms . This 
approach does not give access to a simple analytical expression 
for the perturbation. 

One motivation for this work is that it has become evident 
that the traditional description of the potential perturbation in 
terms of a sine function in the azimuthal direction is not satis- 
factory. It is simple to estimate the stellar density in a region 
where the orbits are crowded and to compare it with the average 
field, as well as to perceive that the perturbations are much bet- 
ter described as narrow grooves or channels in the gravitational 
potential with an approximately Gaussian profile. We present an 
analytical description of the perturbing potential, which is real- 
istic when compared to observations and produces better self- 
consistency compared to other models. Strictly speaking, we 
give priority to the analysis of one aspect of self-consistency, 
which corresponds to the shape of the arms or equivalently to 
the position of the density maxima in the plane of the galaxy and 
to their profiles. We employ the expression "self-consistency of 
arm shape" to describe our results, although several authors that 
we mentioned above simply refer to "self-consistency" for simi- 
lar analyses. Other advantages related to the expression proposed 
in this work for the perturbation are 1) avoiding fluctuations be- 
tween positive and negative values for the potential (as produced 
by a sinusoidal function) over the whole disk for the whole range 
of parameters, 2) controlling the arm-interarm contrast and, 3) 
dealing of the thickness of the spiral arms. 

The focus in this work is on normal galaxies in which the 
structure is dominated by spiral arms, not by a bar. Naturally, 
many galaxies present a weak bar in their center, with an ex- 
tended region of spiral arms beyond the end of the bar. The 
present model can be useful to interpret their structure. 

In addition to the motivations mentioned above, there has 
been considerable improvement in recent years in the knowledge 
of the structure of our Galaxy and, in particular, of its rotation 



curve and pattern speed (discussed, respectively, in Sects. 14.11 
and !6.U . Therefore it is worth constructing a model that is ade- 
quate for comparisons with observations and other models. We 
are conscious that a number of authors consider that our Galaxy 
has a strong bar, so that our model would not be applicable, but 
in any case, the comparison with observations may give us clues 
to this question. Other models that we have already mentioned 
(A&L, Pichardo et al, Contopoulos & Grosb0l) use the same 
hypothesis of spiral-dominated structure and are directly com- 
parable. 

The organization of this paper is as follows. In Sect. 12.11 we 
review the classical model for the spiral arms. In Sect. 12. 2J we 
present our new model for the spiral arms discuss the required 
parameters. In Sect. [3] we deduce the relation between the den- 
sity contrast and the perturbation amplitude. We constrain the 
perturbation amplitude to a range of values by means of obser- 
vational evidence for the density contrast. In Sect. |4]and subsec- 
tions, we compute the axisymmetric potential for two different 
rotation curves. First, we use a realistic rotation curve, which 
is obtained by fitting the data of the Galaxy. The second model 
uses a "flat" rotation curve that is modeled by two axisymmetric 
components (bulge + disk) with different bulge-to-disk ratios in 
order to better understand the effect of the bulge component on 
the stellar orbits. In Sect. 15.11 we present the equations of mo- 
tion and the integration scheme, and in Sect. 15.21 we derive the 
density response. In Sect. [6] we discuss the adopted corotation 
radius and the angular velocity of the spiral pattern. The main 
results are presented in this section, discussed separately for dif- 
ferent ranges of radii: between the Inner Lindblad Resonance 
(ILR) and the 4: 1 resonance, beyond the 4: 1 resonance, and the 
connection between the bulge and the inner stellar orbits. In Sect. 
[7] we give our conclusions. 

2. The spiral arms 

2.1. The classical model 

The surface density of a zero-thickness disk can be represented 
mathematically as the sum of an axisymmetric or unperturbed 
surface density Eo(/?), an d a perturbed surface density £[(/?, tp), 
which represents the spiral pattern in a frame that rotates at an- 
gular speed Q p . The azimuthal coordinate at the rotating frame 
of reference is <p = 6 — Q. p t, where 6 is the angle at the inertial 
frame. As usual, the physical surfa ce density is given by the r eal 
part of the following equation (see Binnev & Tremainell2008l) : 

Si {R, 6 - Q p t) = z se >l»«e-n,,t) + MR)] f (2) 

where f m (R) is the shape function, which describes the spiral, 
and E s is a varying function of radius that gives the amplitude of 
the spiral pattern. From Ei, as described by Eq.|2] we derive the 
potential <!></, using Gauss' law, and obtain 



$>i(R,6-Q. p t) = O d e 



i[m(6-Q p i)+f m (R)] 



where 



= - 



27rGE, 



(3) 



(4) 



In the above expression, k is the wavenumber and G the 
Gravitational constant. Using Eqs.[3]and|4] we get 



(5) 



This is a well-known result, found bv iLin et al. I (fl969l) and de- 
rived from the WKB theory (see also Appendix B). 
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Fig. 1. Azimuthal density profile at different radii for m = 2, 
i = 14°, and cr = 4.7 kpc. Here we set S Jo = 1, since we only 
want to show how the density profile varies with the azimuthal 
angle. 



Fig. 2. Map of the perturbing potential in the plane of the galaxy. 
The colors represent the values of <5i (R, <p) in arbitrary units. The 
picture on the left is the potential with cr — 2.5 kpc and on the 
right cr = 4.7 kpc. For i = 14° we have on the left cr ± = 0.6 kpc 
and on the right <x ± = 1.1 kpc 



2.2. The new potential for the spiral arms 

Before we begin, it is important to clarify what we under stand by 
spiral arms. The model we adopt is based on the idea of Kalnais 
d 19721) of stellar orbits of successively increasing radii in the 
disk, organized in such a way that they get close together in some 
regions, thus presenting excesses of stellar density. 

Traditionally, the perturbed surface density has been de- 
scribed by Eq. [2] This approach, however, is not realistic. The 
brightness profiles observed in galactic disks in circles around 
the center are not sine functions, nor are the density profiles 
obtained theoretically from the crowding of stellar orbits. The 
surface density we consider to be realistic is given by a density 
excess that follows a logarithmic spiral, with a Gaussian profile 
in the transversal direction. The potential that produces this kind 
of density has the form 



<&i(R,tp,Z) = -£oRe~ 



r[l-cos(mp-f m (R))]-s s R-\kz\ 



(6) 



where £q is the perturbation amplitude, ej 1 the scale length of the 
spiral, cr the width of the Gaussian profile in the galactocentric 
azimuthal direction, k = Rt ™ (/) the wavenumber, and f m (R) the 
shape function. The scale-length of the disk and of the spiral are 
not necessarily the same. This question is discussed in Sect. [3] 
The shape function is 



/»(*) 



tan(/) 



ln(R/Rj)+y, 



(7) 



where m is the number of arms, i is the pitch angle, Rj is the 
point where the spiral crosses the coordinate x = 0, and y is only 
a phase angle. 

Solving Poisson's equation using Eq.|6] with the assumption 
of zero-thickness disk and of tightly wound spiral arms (TWA) 
in the plane z — we have 



-\[\-cos(nup-f m (R))] 



(8) 



where S So is associated with the perturbation amplitude (see 
AppendixlBlfor more details) by the equation 



Z S0 (R) 



(pm R 
2nG cr 2 \ tan z'l 



(9) 



The azimuthal profile corresponding to Eq. |8]is illustrated in 
Fig.Q] We emphasize that cr is the half-width of the spiral arms 
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Fig. 3. Radial force due to the spiral arm in the direction ip = 0. 
The circles in red represent the radial force due to our potential, 
while the crosses in green represent the radial force associated 
with the classical sine function perturbation. Here, positive val- 
ues of the force mean that a star would be pushed outward and 
negatives values that it would be pulled inward with respect to 
the galactic center. In this figure, the number of arms is m = 2, 
the pitch angle i = 14", s s = 0.4 kpc 1 , and for our model the 
arm width cr = 4.7 kpc 



in directions of galactocentric circles (see deduction in Appendix 
lAl . Because all these circles cross the logarithmic arms with the 
same angle i, the true width in a direction perpendicular to the 
arms is given by cr ± 



o\l = cr sin(0- 



(10) 



This parameter allows us to reproduce arms of different widths, 
as shown in Fig. [2] 

Figure Opiots the radial force produced by th e spiral arms in 
our m odel (red circles, which are very similar to [Pichardo et al.l 
2003) compared with the classical model (green crosses, derived 
from Eq. as a function of Galactocentric distance R. In both 
cases the perturbation amplitude was set to be equal to one, since 
we just want to see how the force profile varies and are not inter- 
ested in absolute values. The first thing we can notice in Fig. [3] 
is that for the classical model the radial force is stronger for the 
same perturbation amplitude, mainly in inner regions. Thus, the 
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perturbation amplitude in the classical case would be lower than 
in our model. This is due to the fact that a sinusoidal potential 
varies between positive and negative values and that the force, 
proportional to the variations of potential, is stronger than in a 
case where the potential ground level is equal to zero. Another 
difference between the two models occurs around 9 kpc, where, 
in the classical model, the force becomes positive while in our 
model the force slowly tends to zero, but it is still negative. This 
happens because a star in that region feels the effect of the last 
piece of the arm, as we can see in Fig. [2] in the x direction, so 
that it would be pulled inwards. 

3. Relation between density contrast and 
perturbation amplitude 

Following Antoia e t al.l d201 ll) we take as a measure of the den- 
sity contrast 

A 2 ~-^, (11) 

where £</ is the axisymmetric surface density and S so is the max- 
imum density of the spiral arms. This relation is valid under 
the assumption of a mass-to-light rati o of the order of 1 when 
A? is measured in the infrared bands (lKenjll992h . lAntoia et al.1 
d201 lb collected in the literature the estimations of Galactic and 
extragalactic density contrasts, and found them to be in the range 
0.13 < A 2 < 0.23. Therefore, an average value would be about 
A 2 =0.18. 

Table 1. Adopted spiral arms properties 



Property Symbol Value Unit 



Number of arms 


m 


2 




Pitch angle 


i 


14" 




Half width 




4.7 


kpc 


Scale length 




2.5 


kpc 


Spiral pattern speed 




23 


km s" 1 kpc~' 


Perturbation amplitude 




600 


km 2 s~ 2 kpc -1 



The maximum density of the spiral arms is given by Eq. [9] 
We refer to the maximum at a certain radius, not a maximum 
over the whole disk. To simplify the analysis, we assume that the 
axisymmetric surface density has an exponential behavior (Eq. 
[T2b . Galactic disks are often represented by exponential density 
laws and, depending on the value of sj, can explain relatively 
fiat rotation curves. We know that such an exponential law will 
not generate an exact fit to the rotation curve, given by Eq. [TTJ 
but it simplifies the analysis and gives us a good hint about the 
value of the perturbation amplitude 

X d (R) = cr Q e- E " R . (12) 
Therefore the density contrast is 

5f° = 2>i ^ e-fe-^* (13) 

"Ed 2nG cr 2 \ tan /|cr 

The density contrast depends not only on the perturbation 
amplitude £b> but also on cr, which describes the arms width. In 
our model, we can set cr as a constant or a function of radius. If 
we assume that the scale length is the same for both components 
(Ed = s s , which is a reasonable assumption because we do not 



want the amplitude of the spiral to drop too fast outwards) and 
compute the value of <tq using the density of the solar neigh- 
borhood, £</(/? ) = 49 M pc~ 2 (see lBinnev & Tremaind l2008, 
Table 1.1), for Rq = 7.5 kpc and s d l =2.5 kpc we get <x = 984 
M G pc~ 2 . Then, using the values from Table[T]for the number of 
arms and the pitch angle, Eq. [13]reduces to 

|o =310 -4^ (H) 

Since we have used, up to now, parameters of the spiral struc- 
ture that are based on studies of our Galaxy, we shall estimate 
the contrast for a galactocentric distance of about 5 kpc, which 
in our model is about midway between the ILR and corotation, 
as discussed in Sect. [6] Numerically we have R » cr, so that Eq. 
[T4lbecomes 

^«3.10-%. (15) 
A/ 

A plot of the density contrast based on this equation is shown 
in Fig. [4] Comparing the values of A 2 from the literature, as dis- 
cussed above, with this figure, we find a range of amplitudes 
of the perturbation between 400 to 800 km 2 s~ 2 kpc -1 , which is 
compatible with the range of density contrasts of Table 1. Taking 
an average value A 2 = 0.18 gives us £b = 600 km 2 s~ 2 kpc -1 . 



0.36 



0.32 - 




^ 1 ' 1 ' 1 ' 1 ' 1 

100 200 300 400 500 600 700 800 900 1000 
i; (km 2 s" 2 kpc 1 ) 



Fig. 4. Density contrast S , fLj as a function of perturbation am- 
plitude for the Milk Way. 



This is an estimate for R » 5kpc. However as we can 
see from Eq. [T4l th e d ensity contrast increase s with radius. 
iKendall et alJ d201 ll) and Elmegre en et al.l (201 1) showed many 
galaxies with a growth of density contrast with radius. If we 
adopt Ed - e s , the density contrast varies as the square of the 
radius, which seems a little too fast in many cases. This prob- 
lem can be attenuated in two ways: 1) the scale length of the 
spiral and of the disk may be slightly different (sd < e,), which 
will produce an exponential decrease, 2) the arms width may 
increase with radius, producing a density contrast with moder- 
ate increase, closer to the results of Kendall et al. For a specific 
galaxy, it would be possible to check the azimuthal density pro- 
file and its variation with radius and to adjust the parameters to 
match the observations. One merit of our model is that it facili- 
tates such adjustments. 
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4. The axisymmetric component 

The main dynamical information of a spiral galaxy is derived 
from its observed rotation curve. In fact, the rotation curve gives 
us the radial gradient of the potential, since we can compute 
<&o(R) by integrating Eq.[T6l 



F - - 



vL(R) 



(16) 



dR R 
where V ro , is the rotation curve. 

4.1. The rotation curve of the Galaxy 

To contribute to a deep understanding of the spiral structure of 
any given galaxy, we must make the effort of introducing the real 
parameters, such as the dimension of the components, the rota- 
tion curve, and the pattern rotation speed. Although there may be 
controversies concerning some of these parameters, we chose to 
illustrate our model of spiral perturbation by applying it to our 
Galaxy. However, we also used a simplified analytical expres- 
sion for the rotation curve, introduced in the next subsection, 
with the aim of reaching qualitatively a general understanding 
of the effect of the the mass of the bulge on the inner orbits. 

The "real" rotation curve of the Galaxy that we have adopted 
(Fig. [3j is quite flat, except for a local minimum at about 8.8 
kpc. The existence of this dip in the curve is revealed, for in- 
stanc e, by the study of the epicycle frequency in the galactic 
disk ( Lepine et al. 2008). This is confirmed in a recent work by 
Sofu e et al.l (l2009l) . which makes use of different data sets, in- 
cluding data on masers with precise distances measured by par- 
allax, to trace the curve. The rotation curves presented in these 
two papers are quite similar, except for minor scaling factors 
due to the adopted galactic parameters, the solar distance from 
the Galactic centers Rq, and the solar velocity Vo. The curve is 
conveniently fitted by two exponentials and a Gaussian 



V rot {R) = 398 e - s/26 - 014/s + 257 < e -*l&W* 

-20e~ [(R ~ 8 - 8)/a8]2 . 



(17) 



300 



200 




R (kpc) 

Fig. 5. Rotation curve of the Galaxy. The red line is the rotation 
curve fitted using Clemens' data (black dots), and the blue line 
is the pattern speed Q. p x R corresponding to an angular velocity 
£l p = 23 km s _1 kpc -1 . The corotation radius is at 8.4 kpc. We 
adopted Rq = 7.5 kpc and Vo = 210 km s . 




2 4 6 



10 12 14 16 18 20 

R (kpc) 



Fig. 6. Axisymmetrical potential as a function of the galactocen- 
tric radius R. 



The first two components contain terms in (1 /R) and (l/R) 2 
inside the exponential function, which produce a decrease of 
their contribution towards small radii. The last (Gaussian) term 
represents the local dip. The interpretation of a similar curve, ex- 
cept for the dip, in term s of components of the Galax y is given by 
iLepine & Lerovl J2000) and is different from that of Sofue et all 
(120091) . The pattern speed is discussed in Sect. 16.11 

Figure [6] shows the axisymmetric potential that results from 
the integration of Eq.[l6] We used the trapezium rule with adap- 
tive step to solve the integral, and we fixed the potential as equal 
to (Do(lOO) = at R = 100 kpc in order to set the arbitrary con- 
stant. 



4.2. A model of a flat rotation curve with a bulge 

Many galaxies have a flat rotation curve, but some of them 
present a peak in the inner part. For instance, as seen in Fig. [5] a 
peak for our Galaxy appears at a radius smaller than 2 kpc, with 
a maximum around 300 pc. The nature of the peak is a subject of 
debate. Some authors consider that it is not due to rotation, like 
iBurton&Liszild 19931) . who believe that there is a strong velocity 
component of gas outflow from the central regions of the Galaxy. 
Others consider that an important departure from axisym met- 
ric, like a triaxial bulge according to lGerhard & Vietril(ll986l) . is 
required to explain the peak. To perform a first approximation 
analysis on the effects of the presence of a bulge on the stellar 
orbits in the inner regions of a galaxy, we adopt the hypothesis 
that the bulge is axisymmetric, as it seems to be, at least ap proxi- 
mately, in most spiral galaxies ( M endez-Abreu et al.l20 10). Also 
in axisymmetric bulges, one can find elliptical orbits in their cen- 
tral regions. Then, even a very small perturbation could give rise 
to an oval structure. Since in this part of the work we are not 
interested in solving the specific case of our Galaxy, but want 
to have a more general idea of the effect of a bulge, we do not 
use the observed rotation curve of our Galaxy. Instead, we use 
a simpler one from which the potential is obtained analytically, 
thus avoiding numerical integrations. F or this study we adopt the 
model for the rotation curve given by IContopoulos & Grosb0ll 



5 



T. C. Junqueira et al.: The spiral structure of the Galaxy 



(1986), in which the presence of the peak is modeled by two 
axisymmetric components (bulge + disk): 



V mt (R) = Vmax yj fb^R exp(-e b R) + [1 - exp(-e d R)]. 



(18) 



Here 



and e d 1 



are the scale length for the bulge and disk 
components, respectively. The relative importance of the two 
components is given by the bulge fraction f b . 

We analyzed the inner stellar orbits obtained with rotation 
curves with different values of f b , that is, with different impor- 
tance of the bulge. The results are discussed later, in Sect. [6] 



5. Stellar orbits 

5.1. The integration scheme 

Once we know the two components of the potential, <t>o(/?) and 
<&i(R,tp), we can derive the equations of motion governed by 
them. Figure Q illustrates the angles and angular velocities used 
in the equations. When the Galactic plane is represented, the mo- 
tion of stars and the pattern speed are shown clockwise, which 
corresponds to negative values of angular velocities. 




Fig. 7. Scheme of the Galactic plane. The yellow dot represents 
a star rotating with speed in a frame of reference (dashed line) 
corotating with the spiral arms at the angular velocity Q p . 



The Hamiltonian for such a system is given by 



H =2 



-a p Ji+Oo(R) + Oi(R,ip), 



(19) 



where p, and J\ are the linear and angular momentum per mass 
unit, respectively. It is important to note that J\ is measured with 
respect to the inertial frame. The equations of motion are 



dR 



dp r 

dt 



J\ dO (R) 



R 3 



dR 



d®i(R,f) 
dR ' 



'tt - II _ o 
dt ~ Ri 



dJi 
dt 



dOi(fl,y>) 
dip 



(20) 
(21) 
(22) 
(23) 



The equations above were integrated using the 6 lh order 
impli cit Runge-Kutta-Gauss (RKG) (see San z-Serna & Carvel 
11994 for more details), with a fixed time step dt = 10 _1 Myrs. 
To find the stable periodic orbits, we applied the method of 



Poincare's surface of section (Poincare 11892b . The Poincare sec- 
tion was fixed where the orbits cross the axis <p = 0°. The peri- 
odic orbits trapped around circular orbits are represented by a 
sequence of points in the phase-space (R,p r ), lying on a curve 
called an invariant curve. 

In practice, what we do is to compute the energy of a circular 
orbit using Eq. [19] Thus this equation becomes 



h(R) = ^5 - *Vb + «£><>(*). 



(24) 



This expression has no perturbation terms. The angular mo- 
mentum per mass unit is given by Jo = RV ror , and Oo(/?) is the 
galactic potential. For each radius, we have a value of energy h, 
then we come back to the Eq. [T9l with H — h. Thus, we have 
to find the initial conditions (R, p r , J\ , <p) that will close the orbit 
for a given h. The last variable can be fixed at <p = 0, and the 
angular momentum J\ is found solving Eq.[19]for a given R and 
p r . Therefore, we only have to deal with the pair of variables 
{R, p r ). The pair that gives us a stable periodic orbit is found at 
the center of the island in the Poincare section. 

The families of periodic orbits that we obtained (starting 
from energies equivalent to circular orbits spaced 0.2 kpc apart, 
from 3 up to 8 kpc) are shown in Fig. [8] The inner part from 2 up 
to 3 kpc was avoided here; we discuss this region in Sect. 16.41 

5.2. Response density 

The response density is computed by considering the conserva- 
tion of flux betwee n two successive orbits in the perturbed and 
unperturbed cases (IContopoulos & Grosb0lHl986l) . In our case 
we consider a series of circular orbits spaced at 0.2 kpc. We can 
imagine in each orbit a large number of stars equally spaced in 
the unperturbed case. The area S of an angular sector, which 
is between two neighboring orbits and two successive stars in 
each orbit, is transformed into S ' when the perturbation is intro- 
duced. The two areas contain the same mass, so that the density 
is proportional to 1/5". There are two reasons for the change 
in distance between neighboring stars. In the radial direction, it 
is because the orbits approach or recede one from the other. In 
the azimuthal direction, it is because the stars travel at differ- 
ent angular velocities along the orbit. Due to the conservation 
of angular momentum, the stars move more slowly when they 
are at larger galactic radius, and get closer to one another, since 
the flux is conserved. This is easy to take into account, since the 
spacing between stars is A<p oc ipr 1 . The response density can be 
written as 



2, 



esp 



T, c 2nR c AR c At 
TRARAip ' 



(25) 



where E c is the density, R c is the radius, and AR C is the spacing 
between two successive orbits, all of them in the unperturbed 
case, while the same quantities in the denominator refer to the 
perturbed situation. The quantity Af is the time spent by a star to 
move from (p to 9 + Atp, while T is the period. In practice, what 
we do is to divide the circles in N sectors, thus making it possible 
to find the positions R and the time for each angle, defined by iA6 
(i = 1..N). Therefore we can find the time At spent by the star to 
move between two adjacent sectors. 

Since we are interested in the position of maximum response 
density for each orbit, Eq. |35]reduces to 



Af 



RARA6 



(26) 
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As the values of S c , R c , AR C , and T do not change in a given 
circle, they can be ignored in Eq. [25] Thus the position of max- 
imum density is the sector where Eq. [26] reaches its maximum 
value. 

We are interested in the position of maximum response den- 
sity for each orbit. The position of the density maxima is shown 
in Fig. E 
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Fig. 8. Here we show a series of closed orbits in the plane of the 
Galaxy obtained using the rotation curve of Eq. [17] The spirals 
indicated in green represent the imposed perturbation; they co- 
incide with the response density (red dots) up to 4: 1 resonance. 
The red circle is the corotation radius (R cor = 8.4 kpc). The per- 
turbation amplitude is - 600 km 2 s~ 2 kpc -1 . The yellow dot 
shows the Sun position. 



6. Results 

6.1. A model for our Galaxy 

Throughout this work, we have adopted the parameters that ap- 
peared to us to correspond to our Galaxy, like those presented in 
Table 1, and the rotation curve. In the next sections, we present 
both general results, which are valid in principle for other spiral 
galaxies, and specific results, which depend on the exact choice 
of parameters. One of the most important choices when we apply 
a model to a galaxy is that of the pattern speed, or equivalently, 
if the rotation curve is known, of the corotation radius. 

In the last few decades, many papers have pointed to a 
corotation r adius of the Galaxy situated close to the orbit 
of the Sun dMarochnik et al.lll972l ICreze & Mennessierf [19731 
iMishurov & Zeninalll999l among others). An additional argu- 
ment is that the classical theory of spiral arms tells us that 
the arms exist between the Inner and Outer Lindblad reso- 
nances and that the corotation falls roughly midway between 
these resonances. The data from the literature collected by 
IScarano & Lepind ([2012) supports the theory, showing that the 
corotation occurs about midway between extremities of the re- 
gion where the arms are seen. In our G alaxy, the spir al arms 
extend from about 3 kpc to 13 kpc (e.g.. iRusseill 120031) . which 
would situate the corotation at about 8 kpc. Recently, a few di- 
rect methods allowed the decrease of the uncertainty on the coro- 
tation radius to a few hundred parsecs. One example is the de- 
termination of the pattern speed by integrating the orbits of open 




Fig. 9. This figure shows a plot of the observed arm structure of 
the Galaxy in the solar neighbourhood. The Sun is at the center 
of the figure, and the distances on the two axes are relative to 
the Sun. The Galactic center is downward and outside the figure. 
Different kinds of objects have been added to better trace the 
structure, all of them shown with a same symbol (gray dots): 
Cepheids, open clusters, CS sources (see text for details). Curve 
a, in brown, represents a stellar orbit at the 4: 1 resonance, and 
curve b an orbit at the 6: 1 resonance. Curve c is not identified in 
terms of resonances, since it is a linear fit to the observed points. 
The other lines (two in red, one in green) are logarithmic spirals 
fitted to the data. 



cluster s back to their birthplaces. As shown by iDias & L enine 
d2005l) . the result, R CR = (1.06 + 0.08)R sun , does not have any 
strong dependence on the adopted rotation curve. A second di- 
rect observ ation is the position of the ring-shaped HI void at 
corotation (Amores et al . 2009), and a third one the position of 
the square-shaped sp iral arm associated with the 4: 1 resonance 
dLepine et alJl201 la ). The existence of a gap in the distribution 
of open clusters and of Cepheids, along with a step in the metal- 
licity distribution, is naturally explain ed in terms of the coro- 
tation resonance (Lepin e"et al.l 1201 Tbh . What we call a direct 
method is one that does not involve any complex model with un- 
certain parameters and therefore generates robust results . This is 
not the case of n-body simulations, or chemical evolution mod- 
els. Based on the results above, we adopt Rcr = 8.4 ± 0.2 kpc; 
accordingly, with Eq.[17] we find Q. p = 23 km s _1 kpc 1 . 

In the present work, we have not investigated the spiral struc- 
ture beyond corotation, a task left for future work. As we discuss 
in the next sections, there is a change in the nature of the arms at 
the 4:1 resonance. Only between the ILR and that resonance are 
the arms explained by the crowding of stellar orbits. 

There are models of the spiral structure that are in conflict 
with the corotation radiu s that w e advo cate. For instance, model 
of the inner Galaxy bv iDehnenl (|2000) suggests that the radius 
of the Outer Lindblad Resonance (OLR) of the Galactic bar lies 
in the vicinity of the Sun. His model does not belong to the cate- 
gory of direct methods, as we defined above, since it is a simula- 
tion and by nature, based on many approximations and uncertain 
parameters. The author does not claim that he presents a precise 
determination of the OLR of the bar and suggests that this radius 
is * 6-9 k pc. A re c ent re view of pattern speeds in the Galaxy is 
given bv lGerhardl (120111) . He recognizes that the pattern speed 
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Fig. 10. Rotation curves and corresponding inner orbits for different strengths of the bulge. At the top are the rotation curves for 
different values of /j (see text). From left to right, the parameter //, is 3.2, 2.2, 1 .2, and 0, respectively. Below are the orbits associated 
with each rotation curve. The blue line is the bulge contribution, the green line is the disk contribution, and the red line is the bulge 
+ disk contribution. 



of the spiral structure is about 25 km s _1 kpc -1 , but favors the 
idea that the bar rotat es at a higher velocity. T he review does not 
mention the work of Ibata & Gilmorel (fl995l) . who determined 
the rotation speed of the bulge, which was found to be the same 
that we adopt for the spir al arms. In a study of the morph ology 
of bars in spiral galaxies, lElmegreen & Elmegreenl (119851) con- 
cluded that in late-type spirals the bar may extend only to the 
ILR. This is the hypothesis that we adopt to compare our model 
with the Galaxy, based on the fact that the solar neighborhood 
contains spiral arms only. 

6.2. The structure between the ILR and the 4:1 resonance 



Many 
formed 

Amaral & Lepine 
2004; Antoja et al 



s tudies of self-consi st ency 
dContopoulos & GrosboH 1 1 9881: 



Pichardo et al. 2003 



have been 
Patsis et"al] 



per- 



_ _ . , M artos et al. 

201 lh . IContopoulos & Grosb0ll dl986h 



already studied the response density and concluded that strong 
spirals do not exist beyond the 4:1 inner resonance. This result 
was confirmed by the other works and also by the present work. 
In Fig. [8] we show the locii of maximum response density of our 
model. It can be seen that there is a very good self-consistency 
of the shape of the spiral arms in the region inside the 4:1 
resonance. This self-consistency is particularly satisfactory in 
the case of our model, since a potential perturbation with a 
Gaussian profile generates a very similar response density and 
response potential. The fact that two well-behaved logarithmic 
spiral arms extend about the 4:1 resonance is a robust result, 
which was encountered by the various works mentioned above, 
although they use a variety of axisymmetric potentials and 
perturbation potentials. It also corr esponds to an observed 
charac teristic of spiral arms, since Elmegreen & Elmegreen 
(1995) observed that most spiral galaxies have two prominent 
symmetric arms in their inner regions, inside approximately 0.5 
^25, where R25 was defined by them as the radius for which the 
surface brightness is 25 mag arcsec~ 2 . These authors suggest 
that the termination of these arms could coincide with corota- 
tion. This robust structure, with a similar aspect in many spiral 
galaxies, is possibly the one that is able to determine the pattern 
speed. That the pattern speeds are found to have a distribution 



concentrated around 24 km s 1 kpc 1 dScarano & Lepindl20l2h . 
in other words, they are not arbitrary. 

We would like to call attention to an interesting effect. One 
can see that the maximum response density starts to be out of 
phase with respect to the logarithmic spiral perturbation near the 
4:1 resonance. This suggests that the deviation from a perfect 
logarithmic spiral in galaxies can be understood naturally as an 
effect of the 4:1 resonance. A consequence of this is the appear- 
ance of bifurcations. They happen because when the response 
density starts to be out of phase with the imposed density, part of 
the matter remains following the imposed density. This results in 
a weaker piece of arm. In galaxies with a structure similar to that 
of the Milky Way, the appearance of bifurcations is possibly as- 
sociated with the 4: 1 resonance, rather than with the corotation. 
This is not a new result, since it was a lready discovered from 
a theoretical appr oach, for instance, by Patsis et al. (1994]) and 
IPatsis et al](ll997l) . This is also a robust result because it does not 
depend on t he details of a given rotation curv e. Observationally, 
the work of Elmegreen & Elmegreen ( 1995) also states that bi- 
furcations start at the end of the two prominent symmetric arms. 

6.3. Beyond the 4:1 resonance 

In our model, there is a clear change in the nature of the arms 
beyond the 4: 1 resonance. There are no more arms produced by 
the concentration of stellar orbits, but only discrete resonant or- 
bits. How could a resonant orbit behave like a spiral arm? This 
may happen if stars are captured within resonances, similarly to 
asteroids in the S olar system, for w hich the trapping mechanism 
was described by iGoldreichl (11965). For instance, in the asteroid 
belt, the 3:2, 4:3, and 1:1 resonances with Jupiter are populated 
by clumps of asteroids. Similarly, the stars captured within res- 
onances could produce a local increase of density and turn it 
into a kind of arm. If this model is correct, when applied to our 
Galaxy, it would predict that some of the arms seen in the solar 
neighborhood have the form of resonant orbits. 

Direct evidence of an arm with the shape o f a 4: 1 resonant 
stellar orbit was found by iLepine et al.1 (1201 lal) using as tracers 
molecular CS sources, Cepheids with a period longer than six 
days, and open clusters with ages less than 30 Myr. The map 
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obtained in that paper is reproduced in Fig. Here, the same 
symbol has been used for all the kinds of tracers, since we are 
only interested in the general structure. It must be remembered 
that nowadays only in the Solar neighborhood we can have a de- 
tailed description of the spiral structure, since interstellar extinc- 
tion does not allow open clusters and Cepheids to be observed at 
more than a few kpc. 

Figure [9] can be seen as a zoom of the upper part of Fig. [8] 
in the region around the Sun. One can observe the presence of 
a structure with an angle of about 90° (labeled "a" and b rown 
in the figure), which was identified bv lLepine et al.l (1201 lal) with 
the 4:1 resonance. The argument was not only the fit presented 
in this figure, but also much more the points observed along a 
straight portion of about 7 kpc in length, in addition to the fact 
that the resonance is exactly where it is predicted to be, based on 
the known pattern speed. Now, we also present a fit of an orbit 
at the 6:1 resonance to many of the observed points (curve with 
label "b", in blue). This greatly resembles a predicted structure 
in Fig. [8j which passes close to the corner of the arm at the 4:1 
resonance and close to the Sun, in addition to presenting an acute 
corner in the second quadrant of the Galaxy. 

Stable periodic orbits attract around them quasi-periodic or- 
bits, as indicated by the contours of the islands in the Poincare's 
surface of section. It seems that in the region of the disk dis- 
cussed here, quasi-periodic orbits are trapped around the sta- 
ble periodic orbits of the 4:1 and 6:1 families, with very small 
deviations from the periodic orbits themselves. As a result the 
quasi-periodic orbits could be considered as building a "thick 
periodic orbit" structure. These arms, which are composed of 
similar quasi-periodic orbits, would have very small dispersions 
of velocities. 

Finally, one can see in Fig. [8] just below the corner of the 4: 1 
resonance, an orbit that crosses it in two points. This resembles 
the arm-like structure with label C in Fig. |9] We consider that 
this observed structure is probably real, due to the large number 
of tracers that fall on it. This suggests an interesting concept. 
Since resonant orbits can cross each other, it is also possible to 
have arm-like structures which correspond to them crossing each 
other. This contrasts with the classical view of spiral arms, in 
which no such intersections exist. 

Of course, the interpretation we give of some local structures 
may not seem convincing. One possible argument against it is 
that the resonant orbits do not seem to be uniformly populated 
by stars all around them. However, the mechanism of the cap- 
ture of stars by resonances deserves deeper study to determine 
if some regions of the orbits should present a larger population. 
The solar neighborhood is rich in structures that have not been 
explained by any model. One could argue that they are not ex- 
plainable because they are due to stochastic events. The idea that 
we propose here, that some of these structures correspond to res- 
onant orbits, can be verified, since at each point of such an orbit 
we have a particular velocity of the stars that are trapped into 
it. However, this is not the scope of the present paper. We did 
not make any adjustment of the parameters of the model in order 
to fit observed structures. Parameters like Vq (which affects the 
rotation curve), the exact depth and location of the local mini- 
mum in the rotation curve, Q. p , and the width, depth, and pitch 
angle of the imposed perturbation could all have been submitted 
to minor variations. Our experience is that this is a lengthy task, 
with minor variations of the structures produced and some re- 
dundancy (different combinations can give similar results). The 
comparison with the observed data is difficult in part due to its 
poor quality. Besides, some fitted models may look better in one 
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Fig. 11. Poincare diagrams and their respective fam ilies (x\,X2) 
of perturbed periodic orbits (see Contopoulos 1975 , for details). 
In both cases we set the energy at R c = 2.2 kpc, and the pertur- 
bation parameters are the same as Table Q] The only difference 
is the bulge fraction fa; at the top we have = 3.2 and below 
we have ft, = 0. 



region of the galactic plane but worse in another. So, we also 
leave this task for future work. 

Like the previous works on self-consistent spiral structure, 
we did not extend the calculation of stable orbits beyond the 
corotation radius. We expect to see resonant orbits as well, but 
the comparisons with real arms will be even more difficult be- 
cause the distances and shapes of the external arms of our Galaxy 
are poorly known. What is known, from an observational point 
of view, is that the spiral arms show strong departures from loga- 
rithmic spirals in the external regions of galaxies. In most cases, 
the arms look more like a seque nce of straigh t segments that 
roughly follow a spiral path (e.g.. ICherninll 1999b . 

6.4. The bulge and the bar 

About 50% of spiral galaxies have a bar, either strong or weak, 
and there are both kinemat ic and photometric evid ences that 
the Milky Way ha s a bar dBinnev & Tremainal2008l Sect. 6.5). 
IShen etal.1 (1201 Oh argue that the Galactic bulge is a part of the 
disk and suggest that the bulge (in reality, a boxy pseudobulge) 
of the Milky Way is an edge-on, buckled bar that evolved from 
a cold, massive disk. 

In this section, we analyze the connection between the spher- 
ical central bulge and the inner stellar orbits. Figure [TOl shows 
the rotation curves from Eq. [18] The scale lengths for the disk 
and the bulge are fixed at ej 1 = 2 kpc -1 and 1 = 0.4 kpc" 1 , 
Vmax - 220 km s _1 . The parameter fo, defined in Eq. [18] as 
"the bulge fraction", is equivalent to a measure of the "bulge 
strength", since the disk component is kept constant (only the 
bulge is varied), or even to a measure of the bulge density, since 
its scale size is kept constant. We will refer to it as the bulge 
strength. We varied ft from 3.2 down to 0. Below each rotation 
curve in Fig.[l0] we show the family of perturbed periodic orbits 
trapped around unperturbed periodic orbits from R c — 2 up to 4 
kpc; see Sect. l5.1l for more details. For a given energy h, which 
correspond to a radius R c , we may have more then one periodic 
orbit. The family of periodic orbits lies at the center of the island 
in the Poincare diagram. We call it family x\ if the center of the 
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island lies at R < R c , and family X2 if R > R c , (see Contopoulos 
Il975h . 

The parameters for the spiral perturbation are the same as 
given in Table [JJ The orbits in the inner region, almost perpen- 
dicular to the more external orbits, belong in FigJTU] to the fam- 
ily X2 of periodic orbit. They are more prominent where we have 
a strong central bulge (//, = 3.2); as the bulge weakens, the fam- 
ilies of orbits X2 disappear, as we can see in Fig. |TT| In the case 
where fb + 0, we have two families of stable orbits (x\,xi). In 
Fig. QJJ we plot the orbits corresponding to families for which 
we have only small deviations from R c and those with a realistic 
distribution of velocities. As we can see in Fig. [TTJon top, the 
family x\ has a large deviation from R c with a velocity in the or- 
der of p r as -150 km s _1 generating a very eccentric orbit, while 
the family X2 is close to R c with p r as 20 km s _1 ; thus, we expect 
that it would be the favorite family to populate this region. It is 
clear that strong bulges favor the X2 family in the inner regions 
of the Galaxy, while without a bulge we have only the family x\. 

A visual inspection of the orbits associated with different ro- 
tation curves shows elongated stellar orbits in the inner regions 
of the disk for rotation curves with the presence of a peak near 
the center, suggesting the formation of a bar with a half-size ~ 3 
kpc. This model shows the existence of a correlation between the 
formation of a bar and the magnitude of the central bulge. The 
peak in the rotation curve becomes visible when fb is relatively 
large (3.2 or 2.2). A strong bulge generates a prominent velocity 
peak in the rotation curve and modifies the inner stellar orbits. 
In other words, the strength of the bulge has an influence on the 
existence of a bar. This bar appears as a consequence of the same 
potential perturbation that we introduced for the spiral arms, so 
that the bar and the spiral arms have the same rotation velocity. 
The type of bar that we are discussing here corresponds to what 
is understood as a weak bar. In principle, this would seem to be 
more compatible with the model of spiral arms presented in this 
work, in which no effect of a bar is taken into account in the large 
"spiral" region that we described. A bar is weak if it generates 
only a small fraction of the total gravitational field around it. In 
such systems, the stellar orbits can be described by the epicyclic 
theory plus a weak driving force due to the periodic motion of 
the bar. In our Galaxy, the mass associated with the spherical 
bulge is much larger than that associated to the bar, a s deduced 
from infrared luminosity (e.g.. lLepine & Lerovl(l2000 i)). One can 
even expect that at a distance of the order of 1 kpc from the cen- 
ter, outside the bulge but under its strong influence, the potential 
is Keplerian, with a dominant mass at the center. In such a po- 
tential, the orbits are elliptic, with the center of the Galaxy at 
one focus of the ellipse, like the orbits of comets in the solar 
system. In epicycle theory, these are 1:1 orbits. Two families of 
such orbits aligned along a straight line, one on each side of the 
bulge, would be in part responsible for the bar. The alignment of 
such elongated, one-sided orbits with the inner extremity of the 
spiral arms would occur naturally, since stars travel slowly when 
they are at the apogalactic parts of their orbits and this facilitates 
synchronization. Of course, we are not presenting here a com- 
plete model of a weak bar, but only arguing that the existence of 
a weak bar in our Galaxy should not be rejected. If we consider 
that the origin of the spiral structure is related to the interaction 
with an external galaxy, it is more logical to think that the bar is 
induced by the spiral structure, and not vice versa. An interest- 
ing possibility is that since the only imposed perturbation here is 
that of the spiral arms, our results support the idea t hat a pseu- 
dobulge could have evolved from a perturbed disk (She n et al.l 
2010)7 



We must add to the list of similarities between our model 
and the real Galaxy the agreement, in Fig. [8] between the orien- 
tation of the elongated inner orbits, where the spiral arms seem 
to begin, with respect to the Sun. This is very similar with the 
observed angle between the S un and the extremities of the bar 
(see, e.g jLepine et al .112011 at) . 

7. Conclusions 

We have presented a new description of the spiral structure of 
galaxies, based on the interpretation of the arms as regions where 
the stellar orbits of successive radii come close together, produc- 
ing large stellar densities. In other words, the arms are seen as 
grooves in the potential energy distribution. Such an approach is 
not n ew (e.g., ICont opoulos & Grosb0l Il986t lAmaral & Lepine 
1997). The innovation is that the potential energy perturbation 
that we have adopted is itself like a groove that follows a loga- 
rithmic spiral, with a Gaussian profile in the direction perpendic- 
ular to the arms. This represents an important step in the direc- 
tion of self-consistency, since this potential perturbation gener- 
ates, by means of the stellar orbits, spiral arms with a similar pro- 
file. In previous classical models, the potential perturbation was 
represented by a sine function (or a sum of two sine functions, 
if four arms are presented) in the azimuthal direction, but the 
response potential resembled the groove-like one of the present 
work, so that t he self-consistency was poorer. One of the clas- 
sical models dPichardo et al.ll2003h does not use a sine function 
for the imposed potential, but its complexity does not allow an 
easy check of self-consistency. 

A new parameter appears in our description, allowing us to 
control the width of the arms. We found a relation between the 
density contrast between the arm and inter-arm region, and the 
perturbation amplitude, under the assumption of an exponential 
disk. 

We confirm the conclusion of previous works that the 4:1 
inner resonance is a fundamental structure of the disk and that 
similar strong arms do not appear beyond this resonance. This 
result is supported by observations, since it is known that most 
spiral ga laxies have two symmetric prom inent arms in their inner 
regions (Elmegreen & Elmegreen 1995). 

The model is applied to our Galaxy, using a description of the 
axisymmetric potential similar to that derived from the observed 
rotation curve, along with the known pattern rotation speed and 
pitch angle. The range of values for the perturbation amplitude 
compatible with observational evidences of the contrast density 
is 400 - 800 km 2 s -2 kpc -1 . Using a density contrast of about 
20%, an accepted average value, the perturbation amplitude £o> 
is equal to 600 km 2 s~ 2 kpc 1 . It is found that the orbits in the 
region between the 4: 1 resonance and corotation do not reinforce 
spiral arms like in the inner regions, but only features similar in 
shape with the periodic orbits existing in this region. If the model 
can validly explain the spiral arms in the solar neighborhood, 
then the dynamics in that region is determined by the motion of 
stars in quasi-periodic orbits. 

Interestingly, a number of observed structures are very sim- 
ilar to the predicted resonant orbits. The similarities of the 
model with the Galaxy include the orientation of the bar, the 
size and orientation of the 4: 1 resonance orbits as revealed by 
iLepine et al.l (12.0 1 lal) . and possibly a resonant orbit with outward 
peaks or "corners" and inverted curvature between the corners, 
situated close to the Sun. 

The idea that the spiral structure can be self-consistent (a 
certain perturbing potential gives rise to stellar orbits that re- 
produce this perturbation) points towards a long-lived structure. 
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This view is in conflict with some r ecent interpretat ion of the 
spiral arms as transient features (e.g.. ISellwoodll20TTl) . but con- 
sistent with the argument that the spiral structure must be long 
lived base d on the observ ation of the metallicity step at corota- 
tion dLepine etal .11201 lbl) . 

Another interesting result of the present model is that a bar 
naturally appears, without imposing any special condition except 
that the spiral perturbation potential extends a little inwards of 
the ILR. We can see elongated orbits in the central regions of 
the Galaxy, which we identify as being a weak bar. We found a 
correlation between the presence of the elongated central orbits 
and the mass of the bulge. The bulge of our Galaxy is a relatively 
massive one, as revealed by the peak in the rotation curve near 
the center, and consequently can generate a weak bar, according 
to this model. This bar would have the same pattern speed of the 
spiral arms. 

We are conscious that there are many works in the literature, 
that are in conflict with our results. They propose, for instance, 
different corotation radii and different values for the pattern 
speed, the existence of several pattern speeds in our Galaxy, or 
that the arms are transient structures. Another idea is that the bar 
rotates with a different speed relative to the spiral arms. There 
is also a fashionable artistic view of the spiral structure offered 
by a NASA site, often reproduced in scientific works, which 
presents the Galaxy as a perfect logarithmic spiral up to large 
radii. Such a structure could not be a transient one. Still others 
consider that swing amplification underlies the basic physics of 
the spiral arms. To the many contradictory ways of understand- 
ing the Galactic structure, our group of research adds one more 
possibility. We see the arms as relatively long-lived grooves in 
the gravitational potential of the disk, whose shape is primarily 
determined by the stellar orbits and not by shock waves in the in- 
terstellar medium. The spiral structure is self-consistent in terms 
of arm shapes in the range of radius between the inner Lindblad 
resonance and the 4: 1 resonance, so that this region is probably 
the one that guarantees the stability of the whole structure and 
imposes its pattern speed. Outside this region, many observed 
structures have shapes coinciding with predicted orbits at reso- 
nances. Since our model predicts the existence of resonant orbits 
with reverse curvature between corners pointing outward (as in 
our Fig. [8J, it is not surprising to observe structures with similar 
characteristics in our Galaxy. Similarly, it is not too surprising 
to find arm-like structures crossing each other. The expected re- 
sponse of the gas of the disk is to get trapped into the grooves 
and to flow along them. 
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Appendix A: Deduction of the Gaussian potential 
following a logarithmic spiral 

A Gaussian in two dimensions with equal standard deviation cr = 
<t x = cr y is given by 

-(-'-'-0> 2 -<v-vo> 2 

f(x,y)=Ae . (A.l) 

In polar coordinates we have 

( x - Xo ) 2 = R 2 cos 2 (0) - 2R R cos(0) cos(or) + Rl cos 2 (0), (A.2) 
and 

(y - yo f = R 2 sin 2 (f?) - 2RR sin(6>) sin(a) + B% sin 2 (6»), (A.3) 



where Rq = ^jx 2 + y 2 gives the position of the logarithmic spiral. 

If we fix the radii and just look at the variation in the azimuthal 
direction, Eqs. IA.2l and lA.3l become 

(x - x () ) 2 = R 2 cos 2 (0) - 2R 2 cos(0) cos(a) + R 2 cos 2 (6»), (A.4) 
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(y - yo ) 2 = R 2 s in 2 ((9) - 2R 2 sin(0) sin(or) + R 2 sin 2 (6»). (A.5) 



and 



Adding the two equations and using the law of cosines we 
have 



(x - x ) 2 +(y- y f = 2R 2 (1 - cos(0 - a)). 



Taking Eq. lA.6l and putting into Eq. IA.ll this becomes 



f(R,0)=Ae: 



- [l-cos(e-ar)] 



(A.6) 



(A.7) 



Now we can set the angle or as a logarithmic function, which 
defines the spiral arm: 



ln(R/Ri) 
tan(i) 



(A.8) 



but Eq. IA.7I describes only one arm. However, it is simple to 
generalize for m - arms. We simply multiply the term (6 - a) 
inside the cosine, in Eq. IA.7l by m: 



f m (R,e)=Ae^ {l - coslm(e - a)]} . 

This is the final equation to describe m - arms. 

Appendix B: Thin-disk approach for the new 
potential 

The Poisson's equation for a thin disk is 
V 2 (Dj = 4ttGS5(z), 

where we can write the potential as follows: 
Oi(R,9,z) = ®i (R,8)e- M . 

Therefore the Poisson's equation has the form 

r q2 -\kz\ 

W = ^ (eHfe ' V ^ 10 + Ol °^ X 



(A.9) 



(B.l) 



(B.2) 



(B.3) 



where V 2 „ is the Laplacian in polar coordinates. Integrating both 
sides of this equation from z = — °o to z — oo, the second term 
on the right-hand side of this equation disappears, then Eq. 
becomes 



" 2 = tf 



" ' de de 2 



d 2 iff 
M 2 

86 ~~ 
and 

d 2 iff 



The derivatives of iff are 

2R l. m 
a"- \ 2 tan/ 

2 



cty _ 2Rl _ 
dR ~ 



1 



= ^2 + 



m~ 



mR 2 

—r sin *, 

cr 1 



3 m 



-2cos^-- ;sin^ 

tan z i I tarn 



2 D 2 



m R 



■ cosx- 



(B.9) 

(B.10) 
(B.ll) 
(B.12) 

(B.13) 



For a fixed R the maximum density occurs when cos^f = 1, 
sin^ = 0. Assuming a small pitch angle i, the most important 
terms contains tan 2 i in the denominator. Thus Eq. IB.6l is reduced 
to 



v i,e*io 



forfR R 



(B.14) 



Thus the maximum density computed from Eq. IB.4l is 



ft 



m z R 



-E,R 



2nGk cr 2 tan 2 i 



(B.15) 



from the classical spiral density wave theory dLin & Shulll964h . 
The spiral arms were modeled with a cosine function where the 
wavenumber k appears, and it is equal to k — ^ ™ n ^ for a logarith- 
mic description. We know that our expression is not a pure co- 
sine law anymore. However, we can locally think of the Gaussian 
density as a wave for which the amplitude tends to zero far away, 
so that we can keep the classical definition of k and use it on the 
equation above. Taking m — 2 we have 



For the potential <t>io equal to 

<£io = -£ Re$ ll - cosx] - s * R , 

with^ = m(9 - a), and a given by Eq. IA.8l 
Solving V^gOio we have that 

where 
R 2 

iff = —r (1 - cos^) + e s R, 



d>i = 3^-+R 
Y dR 



d 2 iff Idtff 



(B.4) 
(B.5) 

(B.6) 

(B.7) 
(B.8) 



y - & R 

nG <t z \ tan / 

If we chose cr - 
iLin etal.l(ll969l) 



a s,R 



(B.16) 



1 



nG I tan /I 



R, we recover the classical formula from 



(B.17) 



Now the density contrast not only depends on the perturba- 
tion amplitude fo> but also has a dependence with cr, which de- 
scribes the width of the arms. In our model, we can set cr as a 
constant or a function of R. 
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